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I. INTRODUCTION 

Macroscopic quantum effects in two-dimensional 
Josephson junction arrays (JJA's) have been ex- 
tensively studied both theoreticall y 1 ' 2 d dJii&JL and 
experimentall y 8 ^ 10 during the last years. The quan- 
tum nature of the phase of a superconducting order 
parameter is reflected in phase transitions in JJA's. In 
nondissipative JJA's the two main energy scales are set 
by the Josephson coupling Ej between superconducting 
islands and the electrostatic energy Ec arising from local 
deviations from charge neutrality. The ratio Ec/Ej 
determines the relevance of the quantum fluctuations 
and when it increases above a critical value, the phase 
order is destroyed and the array turns into the insulator. 
For large capacitive coupling Ec S> Ej the system can 
be modeled by a renormalized classical two-dimensional 
(2D) XY model. In the opposite limit, the energy cost 
for transferring charges between neighboring islands in 
the array is so high that charges tend to be localized. 
While the nature of the classical 2D XY model is 
well understood, its quantum generalization still poses 
unsettled issues. 

Modern fabrication techniques allow one to make ar- 
rays of ultrasmall superconducting islands separated by 
insulators. In such systems the important factor which 
has a profound impact on the ground state of the JJA's 
is dissipation caused by Ohmic resistors shunting the 
junction o 11 ' 12 ' 13 or quasiparticle tunneling through the 
junctions . 14 i 15 i 16 Despite several experiments with 2D 
JJA' a 17 i 18 and superconducting granular filmsi^ exis- 
tence of the dissipation driven transition and critical 
value of the normal state conductance is at least ques- 
tionable. 



Phase diagrams in quantum JJA's with both mecha- 
nisms of dissipation Ohmic and quasiparticle were stud- 
ied theoretically by Zaikin.— Calculations done within 
the framework of the instanton technique reveal a zero 
temperature phase diagram with two dissipative phase 
transitions. The author claims there are regions on the 
phase diagram where disordered phase and the classi- 
cal Josephson effect could take place. Cuccoli et al£L 
presented an analytical study based on the effective po- 
tential approach. They proposed a model in which two 
different relaxation times lead to the conductance matrix 
with resistive shunts to the ground and among islands. 
Despite of these accurate analytical studies the problem 
of a theoretical explanation of phase diagrams in JJA's 
in dissipative environment is still open. 

Our previous theoretical work^ 2 . in which the attention 
was focused on local dissipation effects, predicted the ex- 
istence of the critical value of the dissipation parameter 
a = 2 independent on geometry of a lattice and mag- 
netic field. Other theoretical studie a n i 12 i 23 i 24 suggest a 
rather broad range of the critical values of the dissipation 
parameter a — 0.5,0.84,1,2 which depends on the di- 
mension of the system and mechanism of the dissipation. 
It seems that an unambiguous in experimental measure- 
ment of the critical value of the normal state conductance 
is elusive. Several group o 17 ' 18 ! 25 ' 26 ' 27 ' 28 ' 29 using different 
experimental techniques obtained various critical values 
of a = 0.5,0.8,1. To explain these theoretical and ex- 
perimental difficulties we propose a model in which local 
(caused by shunt resistors connecting the islands to a 
ground) and nonlocal (shunt resistors in parallel to the 
junctions) dissipation effects are considered. 

The purpose of this paper is to investigate phase tran- 
sitions at zero temperature in two-dimensional capaci- 
tively coupled superconducting arrays with emphasis on 
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the competition of local ao and nonlocal ai dissipation 
effects. The detailed phase boundary crucially depends 
on the ratio of mutual to self-capacitances Ci/Cq and 
specific planar geometry of the array^ Aware of that fact 
we consider capacitive matrix CV,- within the range of 
parameters Ci/Cq which can be adjusted to the exper- 
imental samples. We analyze phase diagrams for three 
different lattices: square (□), triangular (A) and honey- 
comb (H) . We want to emphasize that our approximation 
cannot be used for analysis of the Berezinski-Kosterlitz- 
Thouless transitions since it is appropriate only for phys- 
ical systems where long-range order appears. 

The outline of the rest of the paper is the following: In 
Sec. II we define the model Hamiltonian, followed by its 
path integral formulation in terms of the dimensionality 
dependent nonmean-field like approach. In Sec III we 
present the zero-temperature phase diagram results for 
different JJA's geometries. Finally, in Sec. IV we discuss 
our results and their relevance to other theoretical and 
experimental works. 



II. MODEL 

We consider a two-dimensional Josephson junction ar- 
ray with lattice sites i, characterized by superconducting 
phase <pi in dissipative environment. The corresponding 
Euclidean action reads: 



S = Sc + Sj + <5d , 



(1) 



where 
Sc = 



8e 2 



dr I " y V dr 



Ci- 



i,3 

5j = V / drJi-j {1 - cos [fa (t) - <pj (r)]} , 

(id) Jo 

5d = \ E / dTdT ' a ^ ( t - r ') to M - & ( r ')i 2 ( 2 ) 

id 

and r is the Matsubara's imaginary time (0 < r < 
l/fceT = T is temperature and fee the Boltzmann 
constant (H = 1). The first part of the action §2$ de- 
fines the electrostatic energy where Cy is the capaci- 
tance matrix which is a geometric property of the array. 
This matrix is usually approximated as a diagonal (self- 
capacitance Co) and a mutual one C\ between nearest 
neighbors. We can write a general expression for the dj 
in the following form: 



Cij — 



Co + zC\ for % = j 
-Ci for nearest neighbors 



(3) 



which holds for periodic structures in any dimension; z 
is coordination number of the network. The second term 



is the Josephson energy Ej (Jy = Ej for 



J\ 



\d\ 



and zero otherwise) . The vector d forms a set of z lattice 
translation vectors, connecting a given site to its near- 
est neighbors. The Fourier transformed wave- vector de- 
pendent Josephson couplings Jk are different for various 
lattices. The third part of the action 6>d describes the 
dissipation effects and (r — r') is a dissipation ma- 
trix. We choose two independent damping mechanisms, 
the on-site and the nearest-neighbor, because usually, the 
damping is described in terms of shunt resistors i?o con- 
necting the islands to a ground and shunt resistors in 
parallel to the junctions related to R\. We can write 
dissipation matrix similar to Eq. J3]) in a more closed 
form: 



Oii 



(ao + za\) Si 



a\ y y 8, 



d+d 



(4) 



with the vector d running over nearest neighboring is- 
lands. The dimensionless parameters 



ao 



Rq 

Rq 



Oil 



Ri 



(5) 



describe strength of the local and nonlocal dissipation 
respectively, where Rq = l/4e 2 is quantum resistance. 



A. Method 

Most of existing analytical works on quantum JJA's 
have employed different kinds of mean-field-like approxi- 
mations which are not reliable for treatment spatial and 
temporal quantum phase fluctuations. The model in Eq. 
[2] encodes the phase fluctuation algebra given by Eu- 
clidean group E2 defined by commutation relations be- 
tween particle Li and phase Pj operators, 



Pi = 



— „i<t>j 



[Li,Pi] = 



U,P) 
[Pi,Pj] 



-PS 

P^S ■ 

r i °n i 



0, 



(6) 



with the conserved quantity (invariant of the E 2 algebra) 



PP = Pi 



Pi 



1. 



(7) 



The proper theoretical treatment of the quantum JJA's 
must maintain the constraint in Eq. \7\ A formulation 
of the problem in terms of the spherical model initiated 
by Kopec and Jose^i leads us to introduce the auxiliary 
complex field ipi which replaces the original operator Pj. 
Furthermore, relaxing the original "rigid" constraint and 
imposing the weaker spherical condition: 



N. 



(8) 



where N is the number of lattice sites, allows us to im- 
plementation the spherical constraint: 
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[V4>] 



,-Sc+-d[4>] 



xS[lmipi - P yi (<£)] . 



(9) 



where [Dip] = ^V^Vip* and [V(f>] = Ui^i- I* is con- 
venient to employ the functional Fourier representation 
of the S functional to enforce the spherical constraint in 
Eq. ©: 



S[x{r)] = J 







1 —ioo 





e tf dr\(r)x(r) 



(10) 



which introduces the Lagrange multiplier A (r) thus 
adding a quadratic term (in ip field) to the action in Eq. 

. The evaluation of the effective action in terms of the 
tp to second order in tpi gives the partition function of 
the quantum spherical model (QSM) 



2qsm- / P^]S 



where the effective action reads: 



N e 



(11) 



rP 

Sty] = V / dTdT'HJiji^Sir-T') 

+ W^ 1 (r, t') ~ A (t) SijS (r - r')] 

+ NX(t)S(t-t')}. (12) 

Furthermore, 



Wij(T,T') = / [p0] e l hM-^( T ')] e - 5 c+D[0] ; ( 13 ) 
Zq J 

is the phase-phase correlation function with statistical 
sum 



Z = J [P0]e~ 5c+DM , 



(14) 



where action 6>c+d [4>] is just a sum of electrostatic and 
dissipative terms in Eq. J2j. After introducing the 
Fourier transform of the field 

+ oo 

^( r ) = ^E E v'^" 1 " (is) 

k n=— oo 

with u; n = 2nn/P, (n — 0, ±1, ±2, ...) being the Bose 
Matsubara frequencies. From Eq. lfl3|) the phase-phase 
correlation function reads: 



i 1 v-^ 1 - cos fw„ (t - t') 
>V (r, r) = exp N 



(16) 



The charging energy parameter entering Eq. (fl6|) is 



K 



^(C +4Ci) V^o + 4C : 



4Ci 



(17) 

where K (x) is the elliptic integral of the first kind^ 2 -. 
Furthermore we introduce quantities Eq = e 2 /2Co and 
Ei = e 2 /2Ci related to the island and junction capaci- 
tances. 

The dissipative parameter a and may be explicitly 
written as 



1 



n—>oo N «o + za>i — 2ai£k 

k 



(18) 



where £k is a dispersion and has different form for various 
lattices. In the present paper we consider three different 
geometries of the lattice: square (□), triangular (A) and 
honeycomb (H): 

£^ = cos k x + cos k y , 

££ = co$k x + 2 cos cos ( ~W~k v J , 



with the lattice spacing set to 1 . The results of the sum 
over wave vectors in Eq. [18] are placed in Appendix B. 
Finally, for small frequencies, ao < 2 and cti < 1 the 
inverse of correlation function ([TBI becomes: 



W-UuOn)^ { 8E?" ; n + &&K| fOT W «^ 

71 \ otherwise 

(20) 

In order to determine the Lagrange multiplayer A we ob- 
serve that in the thermodynamic limit (N — > oo) the 
steepest descent method becomes exact. The condition 
the integrand in Eq. (fTTjl has a saddle point A (r) = Ao, 
leads to an implicit equation for Ao: 



(21) 



k n^O 



where 



G 1 (k,w n ) = A - Jk 



8Ec 



2irEj 



(22) 



The emergence of the critical point in the model is sig- 
naled by the condition 



G~ l (k = 0,w n = 0) = A - J o = 



(23) 



which fixes the saddle point of the Lagrange multiplier 
Aq within the ordered phase Aq = Jo- 
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Figure 1: Zero-temperature phase diagram for the total charg- 
ing energy Ej/Ec vs parameter of dissipation an (i = 1 if 
a = and i = if «i = 0) for triangular (A; af tt = 2/3), 
square (□; otf lt = 1) and honeycomb (H; = 4/3) lat- 

tice. Insulating (superconducting) state is below (above) the 
curves. 



III. PHASE DIAGRAMS 

A Fourier transform of the Green function in Eq. i|22l 
enables one to write the spherical constraint (|2~Tj) explic 
itly as: 
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Figure 2: Zero-temperature phase diagram for square 2D 
JJA's with self Co = e 2 /2E and mutual Ci = e 2 /2£i ca- 
pacitance (Eq. I17P for two values of local and nonlocal dis- 
sipation parameter q?o = a i = and Q?o = 2 and ai = 1 
(see Appendix). We can distinguish three areas: phase coher- 
ent state (PCS) where phases in the islands are well defined. 
Insulating state (IS) which could be driven to the phase co- 
herent state by effects of the dissipation (PCSdiss)- Finally, 
insulating state, where superconducting phase is perturbed by 
strong zero point quantum fluctuations due to Coulomb block- 
ade that localizes charge carries to the islands. However sys- 
tem can be driven to the phase coherent state (PCSdils) but 
only by critical values of the dissipation parameters (off 1 — 2 



and a\ r 



!)• 



By substituting the value of Ao = J m ax where J max 
denotes maximum value of the spectrum Jk, and after 
performing the summation over Matsubara frequencies, 
in T — > limit we obtain the following result: 



P. 



f + OO 



p(0 



8E C 



where 



E., 



(24) 



(25) 



is the density of states. We can easily see that solu- 
tion of the model requires the knowledge of the DOS 
for a specific lattice with the superimposition of the self- 
consistency condition for the critical line in Eq. ((24j) . 
A Josephson-junction array network is characterized for 
different lattices by the nearest-neighbor Josephson cou- 
pling Ej with the following wave-vector dependence 



l u - 
7 A - 

J\r — 



Ej£^ 



(26) 



where S^s are given by Eq. [TjJl The Fourier transform 
of the capacitance (dissipative) matrices for a triangular 
and honeycomb lattice^ can be also found in Appendix 
C. 



1 = 



x In 
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2tt 



2E C 



(27) 



The critical values of the nonlocal dissipation param- 
eters have a source in low temperature properties of 
the JJA's correlation function in dissipative environment 
(Appendix A). The dependence of the critical value a.\ 
depicted in Fig. Q] is a direct result of the divergence 
this phase-phase correlator. The Fig. [2] and Fig. [3] point 
out the big difference in values of the self Co and mutual 
Ci capacitance and competition between various dissi- 
pation mechanisms have a severe impact on phase dia- 
grams. In typical real situations mutual capacitance can 
be at least two orders of magnitude larger than the self- 
capacitance what indicates the samples are placed very 
close to Ex/Ej axis in Fig. El 

JJA's devoid of dissipation effects can be in two phases: 
insulator phase (IS) and phase coherent state (PCS). 
However coupling system to the environment we are able 
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Figure 3: Zero-temperature phase diagram for square 2D 
JJA's in space of local ao and nonlocal cti dissipation pa- 
rameters for several values of the ratio Ci/Cb. From the top 
Ci/Co = 1.25,1.67,2.5,5,10,50. Insulating (superconduct- 
ing) state below (above) the curves. 



to drive arrays into PCS even if localization of the charge 
carriers due to Coulomb blockade is strong and dominates 
properties of the system. Furthermore for each geometry 
of the lattice there are critical values of the dissipation 
parameters that lead arrays to situation (Fig. [2]) where 
phases in the islands are well defined and quantum fluc- 
tuations do not perturb a superconducting phase - region 
describe as PCS^ll- Notwithstanding between these two 
boundary situations there is a region on the phase dia- 
gram in Fig. [2] where concrete situation depends on the 
values of the parameters. In this area PCSdiss system 
can be driven to the PCS but coupling to the environ- 
ment does not have to be so strong as in -PCS^™' case. 



IV. RESULTS 

Until now only three papers considered effects with 
both mechanisms of the dissipation . 11 v 20 , 21 The most in- 
teresting is Cuccoli's work where authors introduced the 
full conductance matrix for triangular and square lat- 
tices. It seems their results improve the quantitative ac- 
curacy; nevertheless problem of the theoretical explana- 
tion of the phase diagram of JJA's in dissipative environ- 
ments is thus open. 

A model of an ordered array of resistively shunted 
Josephson junctions was also considered by Chakravarty 
et. alJ^ and simplified at several points. They assumed 
that capacitance matrix is diagonal Cy = CSij . The au- 
thors claim the results do not depend sensitively on de- 
tailed form of Cij. Moreover the matrix ay = h/4e 2 Rij 
where Rij is the shunting resistance between grains i and 
j is reduced to the form in which the information about 
the geometry of the lattice is not included. The obtained 
zero-temperature phase diagram reveals the fact that the 
critical value of the dissipation exists and is proportional 



to the inverse of the dimension of the system which gives 
us critical value a = 1/2 for a square lattice, but espe- 
cially at low temperatures, variational methods are not 
precise enough to perceive such a subtle transition. 

In our model the critical value of the nonlocal dissipa- 
tion parameter ct\ behaves similarly. It depends on the 
maximum value of the Jk spectrum. Because Jk exhibits 
different characters for various lattices hence we could 
observe phenomenon such as nonmonotonic dependence 
of the critical value of the nonlocal dissipation parameter 
for various geometries of the array (see Fig. Q]). When 
we assume diagonal form ay = ao#y then obviously our 
results will not change when we change the geometry of 
the lattice simply because the shunt resistors connecting 
the islands to ground are the same for each island. On 
the other hand if we take into consideration that ao = 
and only a± is present, the situation changes because now 
values of the dissipative matrix strongly depend on the 
Jk spectrum which indicates various values of the matrix 
a.ij depend on the geometry of the structure. This case 
is present in arrays with shunt resistors in parallel to the 
junctions. 

A standard way to study models with only diagonal 
charging energies Cy = CSij corresponds to a com- 
plete absence of screening by the other islands in the 
array.— Notice the mutual capacitance can be at least 
two orders of magnitude larger than the self-capacitance 9 
C\ ~ 10 2 Co. We propose a more realistic model in which 
both self and mutual capacitances are nonzero. To see 
how significant this consideration is we shall analyze Fig. 
03 If we assume, that C\/Cq ~ 50 (see the lowermost line 
in Fig. [3]) we can see the critical values of the dissipation 
parameters change dramatically from local ao = 2 and 
nonlocal a\ = 1 (Appendix A) obtained for fewer real- 
istic cases Ci/Cq = 1 in which both of the capacitances 
are comparable. In the limit (C\ ^ 0, Co = 0) the lattice 
model is equivalent to the Coulomb gas model with crit- 
ical properties not fully understood at present. However 
we have to emphasize that the range of the Coulomb ma- 
trix becomes infinite when Co is set equal to zero. The 
phase boundary in Fig. [3] shift downward with increasing 
ratio Ci/Cq because lower value of the Coulomb inter- 
action between nearest neighbors reduces strong quan- 
tum phase fluctuations and in consequence we observe a 
growth of long-range phase coherence. 

Yagi et alM- experimented with the superconductor- 
insulator transition in two-dimensional network of 
Josephson junctions in detail by varying the junctions- 
area. It was observed the critical tunneling resistance 
exhibited significant junction area dependence. The 
low-temperature behavior of the total charging energy 
Ej/Ec as a function Rq/R u , where R n is the tunnel- 
ing resistance exhibits the same behavior as the curve 
obtained from our theory for square lattices (see Figure 
[T]) in the absence of the local dissipation effects ao = 0. 
The observed phase boundary, which is bending down- 
ward, and the critical value of the nonlocal dissipation 
parameter a™' = 1 is in excellent accordance with our 
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results. 

In order to investigate the effects of quantum fluctu- 
ations and dissipation in JJA's, another group^ made 
two-dimensional arrays of small junctions with various 
Ej/Ec and resistors which caused dissipation effects. 
The value of Ej was controlled by varying the tunnel 
resistance. Each island was connected to the neighbor- 
ing ones by a shunt resistor as well as the tunnel junc- 
tion. The resistance of the shunt resistors was tuned by 
varying their length. Ground states of 2D Josephson ar- 
ray in Ej /Ec-Rq/Rs parameters space reveal the same 
behavior as previous experimental results but there is a 
difference in critical value of the Rq/R s ~ 0.5 which also 
differs from value a™* obtained in this paper. These dis- 
crepancies between experiments can be explained in the 
framework of our model by taking into account differ- 
ent value of the ratio junction-to-self capacitances which 
reduces the critical value of the Rq / R s and considering 
that not only nonlocal dissipation mechanism is present 
an ^ we can obtain Rq/R s ~ 0.5 value (see Figure [3]). 
Real numbers strongly depend on the properties of the 
junctions used in experiments. 



DOS A H □ 



J max / E J 



A 2 



Table I: Maximum values of the spectrum J (k) for three dif- 
ferent geometries of the lattices: triangular (A), honeycomb 
(H) and square (□) 
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Appendix A: SOME PROPERTIES OF THE 
CORRELATOR 



Assume that ao = we write expression for the phase- 
phase correlation function (similar to equation used in a 
previous calculations^ but modified by dissipative ma- 
trix) in form: 



W(r)=exp{ 1^ 



1 — COS {ljJ n T) 



2tt Ej 



(Al) 



V. SUMMARY 

We have calculated quantum phase diagrams of two- 
dimensional Josephson junction arrays using the spher- 
ical model approximation. The calculations were per- 
formed for systems using experimentally attainable ge- 
ometries for the arrays such us square, triangular and 
honeycomb. The ground state of the Josephson coupled 
array with a triangular lattice appears to be most stable 
against the Coulomb effects. This geometry is also the 
case in which the global coherent state emerged when 
the value of the nonlocal dissipation parameter a% is the 
lowest. In JJA's we can observe the phase coherence 
transition which is caused by electrostatic and dissipa- 
tive effects. The detailed phase diagrams crucially de- 
pend on the ratio junction-to-self capacitances, Ci/Cq 
and both dissipation mechanism have a big impact on 
phase boundaries. The nondiagonal terms in capacitive 
and dissipative matrices can change the phase diagrams 
of the system drastically . It is necessary to take them 
into considerations when we have different sources of dis- 
sipation such us shunt resistors connecting the islands to 
a ground and shunt resistors in parallel to the junctions. 
The experimental observation of an universal resistance 
threshold for the onset of the global coherent state seems 
possible, but appears to be difficult. 



It is easy to see the sum over u> n is symmetric when we 
change u„ — > — w n . The key to obtain the solution is 
a calculation the sum or the integral under the expo- 
nent in Eq. (|Aip . Because we are going to investigate 
low-temperature properties of the correlation function we 
could write jj XL,, ~~ * ~kw I-oo ( ^ L0 - ^ n th & t case (Siting 
rid of abs) for large value r we write 



W(r) = exp 



exp 



1 



dtu- 



1 — cos (tui) 



2jEj_ 



ME" 

aiJkE c 
AttEj 



ai Jk 



2tt Ej 

2Ej/axJ k 

ri I (A2) 



where 7 = 0.57721 is the Euler-Mascheroni constant. 
Finally, after Fourier transform we see that correlator 
W _1 (uj m ) ~ |cj m | 2£ ' J / Ql ' 7mBX ~ 1 at zero temperature di- 
verges for ai > 2_Ej/J max . Quantity J max /Ej means 
the maximum value of the Jk which differs for various 
lattices (see Table [J). 



Appendix B: DISSIPATION PARAMETER FOR 
CONSIDERED LATTICES 



In this appendix we give the explicit formulas for the 
dissipation parameter discussed in Sec. II and III. 
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where 



K(x) 



ir/2 



o y/l — x 2 sin 2 <j) 



(B2) 



is the elliptic integral of the first kind 3 - 2 - and the unit step 
function is defined by: 



Q(x) 



1 for x > 
for x < 



(B3) 



4 1 / 4 (2<) 1/2 
(2t- 1) 3/2 (2i + 3) 1/2 



(Bll) 



with t = (a + 3«i) /2a!. 

Appendix C: DOS FOR CONSIDERED LATTICES 

In this appendix we give the explicit formulas for the 
density of states discussed in Sec. II and III. 



For small values of the a± we can write dissipation pa- 
rameter for square lattice as: 

a n =a + 3ai--— + -^+0(af), B4 

for large values values of the ai: 
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2. Triangular lattice 
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1. Square lattice 



AO 



:K 4 1 



e i 



2. Triangular lattice 



(CI) 



where 



a A = — -= K (kj 
7rV3 ai 



(B6) 
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where 



_J_ K ( 
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7T V K 



© ( £ + 2 ) - « - 3) 



(C2) 



1/9 i 3 / 2 r 1/? i 1 / 2 
(2i + 3) 1/2 -l (2t + 3) 1/2 + 3 



4(2t + 3) 1/4 



(2i + 3) 1/2 - lj [(2i + 3) 1/2 + 3 
with t = (ao + 6ai) /2ai. 



1/2 



(B7) 
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©ie+2 ) - (e + i: 



4^3 + 2^ [e (e + i)-e(e-3)] 



(C3) 



Q(C+2 ) - e ^ + 1 ) 



Ki = 4^3 + 2£ 

'3 + 2^3 + 2^-^) [© (e + l)-©(C-3^C4) 



3. Honeycomb lattice 



a 



_x 1 g ao + 3a 



iK(«) 



where 



(2<-l) 3/2 (2i + 3) 1/2 



(B9) 



(BIO) 



3. Honeycomb lattice 



P H (0=4|C|p A (3-4£ 2 ) 



(C5) 
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